Dynamics of a two-component atomic Bose-Einstein condensate 
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The dynamical population oscillations between two internal states of a 
Bose-Einstein condensate are investigated within the rotating wave approxi- 
mation. Analytical expressions for the population imbalance in the number 
states and coherent states have been derived, which predict different revival 
periods. Thus the true quantum state of the condensates may be unambigu- 
ously determined by detecting the atom intensity evolution for one internal 
state. 
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The experimental observation of the Bose-Einstein condensation of a trapped, dilute gas 
of alkali atoms, and the high accuracy of the engineering, are opening a new avenue to in- 
vestigate the interplay between macroscopics and quantum coherence. In particular, recent 
experiments on two-component Bose-Einstein condensates (BECs) in 87 Rb [|],[| have stim- 
ulated considerable theoretical work on the dynamics of the phase and number fluctuations 
of Bose condensates. Vigorous efforts have concentrated on the issue of temporal phase co- 
herence between two coupled Bose condensates. Along semiclassical (mean) analysis much 
attention has focused on the coherent phenomena such as Josephson effect and macroscopic 
quantum self-trapping (MQST) 0,fflKlH. Studies involving small quantum corrections show 
that due to the nonlinearities arising from the atom-atom interactions and the discrete spec- 
trum of the many-body system, the dynamics of the population oscillations arising in the 



mean-field approximation is modulated by the collapses and revivals [|,^,||,||,|T0|,|TT],|T2|.|T3[ . 

In the present work, we study quantum dynamics of a two-component atomic BEC such 
as 87 Rb. Our purpose is to provide an analysis on how the coherent atom oscillations 
between two components depend on the initial states chosen for the condensates. We chose 
two different initial states of the condensate. In the first example we consider Fock state for 
each component of the condensate. In this case, we show that the revival period depends 
on whether the total atom number N is even or odd. In the other example we invoke Bose- 
broken symmetry and consider the product of two coherent states. In this case we show 
that the revival period is always definite. In experiments, such as 87 Rb, initial Fock states 
can be created by condensing independently two thermal atomic clouds trapped in different 
Zeeman level. Whereas, the initial coherent state may be created by splitting a condensate 
(trapped, say, in the Zeeman state \f = 2,rrif = 2)) in two condensates trapped in the 
Zeeman states \f = 2,mf = 2) and \f = 2,rrif = 1), respectively. Another aspect we are 
interested in this paper is the dependence of the coherent atom oscillations on the relative 
phase between two internal states and nonlinear atomic collisions. 

We consider here Bose-Einstein condensation of a trapped gas of atoms that have two 
internal states |1) and |2) with energies Hu /2 and —Huj /2, respectively. There is a spatially 
uniform radiation field with frequency u e that couples the two internal states with a Rabi 
frequency Qq. The atoms in states |1) and |2) are subject to isotropic harmonic trapping 
potentials V^(r) = Ima^r 2 for i = {1,2}, respectively. Furthermore, the atoms interact via 
elastic two-body collisions through the interaction potentials Vij(r — r') = 47r ^ a »j §( v _ 
Where ay is the s-wave scattering length between atoms in states i and j. It is assumed that 
ciij > corresponding to repulsive interactions. In the formalism of second quantization, 
the system is described by the Hamiltonian 
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Here, the atomic field operators ^(r) and $f(r) have been written in the field interaction 
representation which is rotating at the frequency of the external field cu e , 5 = u — u e denotes 
the field detuning from resonance excitation. 

Now we can approximate the field operators in the zero-temperature limit by a two-mode 
model such that \l/j(r) q0j(r), where q is the annihilation operator which obeys the usual 
boson commutation relations. In the two-mode approximation the Hamiltonian becomes 
(with h = 1), 



where 
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The analysis of Eq. (2) is greatly simplified by the introduction of the angular momentum 
operators 
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J y = -cfc 2 ), (4b) 



Jz = ^(c|c 2 + cjci). (4c) 



The Casimir invariant is therefore 
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where iV = cfdi + c^c.2 is the total number operator, which is a constant of motion and 
thus can be set equal to the total number of atoms N. Therefore the present two-mode 
approximation is transformed into an angular momentum model with total angular momen- 
tum given by j = N/2. In terms of angular momentum operators the Hamiltonian (2) may 
be rewritten as 

H = nj z + 2K.r x + f(j), (6) 

where 

-{— + — - K12), (7a) 

« = 2(^- + ^- + K i2) ? ( 7b ) 

/(j) = 2<? 2 + (£1 + £ 2 - « - k% (7c) 

In deriving Eq. (6) we have assumed that 

£i + <5-£ 2 + (^-^)(2j-l) = 0, (8) 

a condition that can always be achieved by adjusting field detuning S and shifting the energy 
levels of the condensate components. Note that some control of the nonlinear parameter 
k can be achieved through the proper engineering of the trapping potential or tuning the 
scattering lengths via Feshbach resonances. 



In general the Hamiltonian (6) is difficult to treat in an exact way because of the presence 
of nonlinear atomic interactions. In order to gain physical insight into the dynamics of such 
a two-component BEC problem, some approximations are necessary: a common assumption 
]T^,|T5[ is the so-called rotating wave approximation (RWA) U = e~ intJz . Using the operator 
identity 

exp(iA J z ) J x exp(—i\J z ) = J x cos A — J y sin A, (9) 

the Hamiltonian in the rotating frame (RF) becomes 

H RF = U + HU + i^ r U (10) 
= ± K ( e 2int Jl + e-^J 2 + J+J_ + JJ + ) + f(j), 

where J± = J x ± i J y are raising and lowing operators. Now we invoke RWA by suppressing 
rapidly oscillating time-dependent terms in Eq. (10). Thus we obtain the RWA Hamiltonian 
in the rotating frame as follows 

hU = -kJ1 + 9U), (11) 

where g(j) = nj 2 + f(j). Instead of discussing dynamics in the rotating frame, we return to 
the Hamiltonian in the previous field interaction representation 

H eff = UH^ F U + - i^U = -kJI + QJ Z + g{j), (12) 

which is our starting Hamiltonian for the following discussions. It reveals in Eq. (12) that 
J z is now a constant of motion and solely determined by the initial condition. The frozen 
J z behavior can be understood when Q is sufficiently larger than k, i.e., Q ^> k. In this 
case, the external field forces the total spin to remain polarized in the ^-direction because it 
costs energy to change the spin direction. To see the validity of our RWA Hamiltonian, we 
numerically integrate the Schodinger equation with the original Hamiltonian (6) and present 
in Fig. 1(a) the expectation values (J z ) against time for various values of Q with the initial 
angular momentum state \j, —j) x where the subscript denotes the quantization axis along 



the x direction. We remark that for the system starts from \j, —j) x , the only nonvanishing 
spin component is J z because ( J x ) = ( J y ) = at all times. Clearly it reveals in Fig. 1(a) 
that when Q = 0, (J z ) vanishes after some time, which means that the initial coherence and 
polarization of J z is completely destroyed by the phase diffusion due to nonlinear term 2k J%. 
However, when Q is sufficiently larger than k, the value of (J z ) remains almost unchanged 
during time evolution, consistent with Eq. (12) in which J z is a constant of motion. As a 
further illustration we calculate the energy spectrums of the two Hamiltonians and present 
the results in Fig. 1(b). It is found in calculations that the difference between the energies 
of two Hamiltonians approaches to be negligible when increasing the value of Q. Thus we 
expect that the RWA Hamiltonian is valid in the parameter regime Q ^> jk. In the following 
we will neglect the constant term g(j) because this term denotes a simple constant energy 
shift and has no contribution to the dynamics of the system. 

Now the Heisenberg equations of motion for the raising and lowing operators (J± = 
J x ±iJ y ) can be readily solved as follows 

j + (t) = e «(-2«J,+«+n) j +j (13a) 

J_(t) = e lt ^ +K -^J_. (13b) 

Given the initial state | ip(0)) of the system, a quantity tailored to the dynamics is the 
population imbalance between the two components 

N-(t) = 2(^(0) | J x (t) | V(0)) (14) 
= e l(K+n)t (4>(0) | e- 2tKtZ J + | V(0)) + c.c. 

Number state description. We shall now proceed to demonstrate how the population 
imbalance evolves for two different kinds of initial states of the condensates. First we consider 
the case that the condensate is in the Fock (number) states, equally the angular momentum 
states | j,m) x where the value of 2m gives the atom number difference between the two 
BECs. We suppose that the condensates are initially localized in one well, corresponding to 
the state | j, j) x . Substituting this initial state into Eq. (14) gives 
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N_(t) = e l(K+n)t yU + m)(j-m + 1) 
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where the matrix elements d 3 mm ,{6) = (j,m \ exp(—i9J y ) \ j,m f ). Using the identity 
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and recursive relation 
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we obtain 



JV-(t) = ^■e i <" +n "e iSj '"' £ (j + m) 
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= Ar[cos(/tt)] JV - 1 cos(fit). 



(18) 



From Eq. (18) one can see that the population imbalance involves a rapidly oscillating 
part and an envelope function that is responsible for the collapse and revival of the population 
imbalance. The most prominent feature in Eq. (18) is that the parity of total atom number 
N, whether it is odd or even, plays an essential role in determining revival periods. When 
iV is odd, the revival period is ir/ k, and the population imbalance is always larger than zero. 
That is, the population undergoes the macroscopic quantum self-trapping. However, when 
the number of atoms is even, iV — 1 is odd and the condensate will undergo antirevivals 
when t = (2n + 1)-k/k so that N-(t) = —N at these times. In this case, the revival period 
becomes 2-k/k and the time average of N_(t) is zero. Two examples are plotted in Figs. 
2(a)-(b). To show the validity of the analytical results of Eq. (18) which is derived from 
the RWA Hamiltonian, we numerically calculate the Schrodinger equation with the original 
Hamiltonian (6). The results are shown in Fig. 2(c)-(d), corresponding to Fig. 2(a)-(c), 
respectively. Clearly, the RWA formula, equation (18), describes the system's evolution 
very well when compared with the exact numerical results shown in Fig. 2(c)-(d), implying 



the different periods of quantum revivals for different parities of N. We notice that the 
reduction in revival amplitude shown in Fig. 2(c)-(d) is due to the state mixing induced 
by angular momentum operator J x , which is eliminated in the RWA treatment. If the Fock 
states with a fixed number of atoms are appropriate for the description of the condensates, 
then by monitoring the initially unoccupied internal state using off-resonant light scattering, 
this change for the magnitude of the revival period of the population imbalance between the 
two components may be observed. A similar phenomenon has been discussed in Ref. [16] in 
the context of the interference pattern, without considering coupling and nonlinear atomic 
collision between the two components. However, because the visibility of the interference 
pattern is never negative, so the revival period considered in Ref. [16] is always n/n with a 
7r phase shift for different parities of N. 

Coherent state description. We turn now to the coherent state description of the problem. 
It is customary to consider the condensates to be in two coherent states, associated with 
a macroscopic wave function with both magnitude and phase, the presence of which is due 
to the spontaneous breaking of gauge symmetry. This is the view we will adopt here, to 
assume the condensate to be such coherent states with definite relative phase 

|^(0)) =!«,/?>■ (19) 



The amplitude of the two components may be presented by |a| = V N sin(9) and \(3\ = 
\J~N cos(0), conserving total atom number. The relative phase between the two com- 
ponents is defined to be 0. Thus the two complex amplitudes are related by f3a* = 
|7Vsin(2#) exp(i(j)). To give a closed expression for Eq. (14) with the initial state Eq. 
(19), we need the normal-ordering arrangement for the operator function exp(AJ^) J + . With 
the aid of the following identity 

—Xcjci \ 1 "_c+r^r 



e - - = £ -cTc[, (20) 

r=0 '• 

where x = e~ x — 1, and after a straightforward calculation, we obtain 

^ 1 oo „m„n,.q~p 
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2=-oo f ' 

x^^X^-cO, (21) 

where y = (e 2A +l)/2, and 2 = (e 2A -l)/2. Substituting Eqs. (20)-(21) into Eq. (14), we have 
the final result for the population imbalance at time t in the coherent state representation 

nt 

N-(t) = Aexp(-2iVsin 2 — ) cos[W sin(/ct) + fit + </?], (22) 

where we have defined A = N^Jl — cos 2 sin 2 (20), W = — iVsin(20) cos0, and a constant 
phase shift </? = arctan[sin0tan(2#)]. We can see from Eq. (22) that the dynamics of the 
population imbalance is mainly determined by the envelope function exp(— 2iVsin 2 y), with 
modulated by a rapidly varying cosine term. The characteristics of collapses and revivals 
can be clearly seen from Eq. (22). For comparison with Eq. (18) derived in the number 
state space, let us chose the value of 9 = 0, corresponding the case that the condensate is 
initially localized in one internal state. Then N_(t) in Eq. (22) reduces to 

nt 

N-(t) = iVexp(-2iVsin 2 — ) cos(ftf)- (23) 

Although no relative phase information occurs in Eq. (23), quite different from number 
states, revivals for the condensate in the coherent states occur at definite times t = 2im/K (n 
is an integer ), independent of the parities of the average particle number N. As illustrated 
in Ref.[16], this difference of revival periods between the Fock state and coherent state 
descriptions originates from the fact that the atom number difference operator c^c-i — cfci 
is quantized in units of 2 when the condensates are in Fock states and units of 1 for coherent 
states. 

We emphasize here that for any realistic experimental set up, the exact atom number is 
unknown and maybe has a large number fluctuation over many different experimental runs. 
As a result, what we obtain in Fock states should be wieghted by a factor p(N) describing 
the probability the state \ j,j) x occurs at the begaining. According to quantum measurement 
theory, however, the atom number in any given experiment would be well defined. Thus 
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in the Fock-state description, the occurence of antirivaval phenomenon, which feaures the 
parity of the atom number, is robust for any single experiment. On the other hand, the 
antirevival phenomenon never occurs in the coherent-state description, which can be seen 
from Eq. (23). This is because the coherent state is a superposition of the Fock states with 
the weighting factor satisfying the Poisson number distrbution. Consequently, in each single 
experiment quantum fluctuations will be present in the coherent state description, which 
smears out the antirevivals. This fact allows us the possibility to detect the true quantum 
state of the two-component BEC system in realistic experiment. 

In the case of the values of nonlinear interactions satisfying the following condition 



our RWA treatment of Hamiltonian (6) is exact and Eq. (22) reduces to description of a 
simple Rabi oscillation with period lir/VL. This form of atomic collective excitation has 
been observed in Ref . [2] . The existence of Rabi oscillation also depends on the initial state 
conditions, which is reflected by the fact that the oscillation amplitude A in Eq. (22) is 
determined by the initial population distribution and relative phase between the two com- 
ponents. In particular when the value of 9 = 7r/4, corresponding to equal atom distribution 
between the two internal states, and = nir, we arrive at the conclusion that the oscillation 
is completely suppressed and no quantum tunneling occurs during time evolution. Note 
that the present conditions of zero oscillation is identical to the fixed points discussed in 
Ref. [5], in which the semiclassical approximation has been adopted. In general the value of 
effective interaction k is not zero and it can be adjusted through the proper engineering of 
the trapping potential, and hence of the condensate wave functions 0«(r). To illustrate the 
interplay between the weak nonlinearity and strong Josephson-like coupling, we present in 
Fig. (3) iV_(t) [Eq. (23)] for several values of k. It reveals in Fig. (3) that the presence of 
the effective nonlinear interaction modulates the fringe visibility of the internal Rabi oscil- 
lations as a manner of collapses and revivals, as has been verified by JILA group 0j. When 
k increases to large values, the rapid phase diffusion caused by the nonlinear interaction 




(24) 



10 



destroys the collective excitation of atoms, corresponding to the complete smearing of the 
Rabi oscillations, as illustrated in Fig. 3(d). 

To see the effect of the initial relative phase on the dynamics of the system, we present 
in Fig. 4 the time evolution of N-(t) for four different values of 0. The value of 9 in these 
figures is chosen to be 7r/4, corresponding to the two condensates with equal atoms. We can 
see in Fig. 4 that the oscillation amplitude is maximal for the values = n/2. Experimental 
observation of this change of oscillation amplitudes may help us determine the relative phase 
between the two condensates. 

We have restricted our analysis to zero temperature. Our study of the effect of initial 
state selection on the consequent dynamics of a two-component condensate should be ex- 
tended to finite temperatures, where dissipative effects associated with the thermal cloud of 
noncondensate atoms must be included. Further study of damping in this system is needed. 

In summary, with the help of rotating wave approximation, we have derived analytical 
expressions for the population imbalance between the two Bose condensates in the number 
states and coherent states. The time evolution of the coherent atom oscillations is shown 
to be quite different between these two state descriptions. For the condensates in number 
states, the revival periods of the oscillations are either tt/k or 27r/K, depending on whether 
the total atom number is odd or even. Whereas, for the condensates in coherent states, 
the revival periods are always 2ti/k. So the true quantum state of the condensates may be 
unambiguously determined by detecting the atom intensity evolution for one trap. 
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Figure captions 

Fig. 1. (a) Time evolution of (J z ) for the value of total atom number N = 100. The 
initial state is \j, —j) x ] (b) energy spectrums of the exact Hamiltonian (closed triangles) and 
the RWA Hamiltonian (closed squares) for the values of A^ = 100, Q = 30k. 

Fig. 2. Analytical results of population imbalance N_(t) for the value of total atom 
number (a) N = 100 and (b) N = 101. The exact numerical results corresponding to (a) 
and (b) are shown in (c) and (d), respectively. Other parameters are n = 0.001f2. 
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Fig. 3. Population imbalance N-(t) for the value of effective nonlinear coupling (a) 
K = 0, (b) k = 0.001ft, (c) re = 0.005ft, and (d) re = 0.025ft. Other parameters are N = 100 
and 6 = 0. 

Fig. 4. Population imbalance N-(t) for the value of relative phase (a) = 0, (b) 
= 7r/6, (c) = 7r/4, and (d) = 7r/2. Other parameters are N = 100 and k = 0.01Q. 
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